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ABSTRACT 

This thesis developed expressions for the aperature fields and inte- 
gral equations for the far field radiation characteristics of the Parallel 
Plate Luneberg Lens Microwave Antenna operating in the TEM or TE^q modes. 
These expressions were programmed on a digital computer to predict the 
far field radiation patterns for several feed systems. Experimentation 
produced far field radiation patterns that were very close to the theor- 
etical patterns for the TE^q mode lens and substantially different for 
the TEM mode lens. 
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I. INTRODUCTION 



Theoretical work by R.K, Luneberg [Ref. 1] on optics in a medium of 
variable index of refraction resulted in the development of the Luneberg 
Lens which has many applications in microwave antennas. Two of the most 
common Luneberg Lenses are of the spherical and cylindrical types, in 
which the index of refractions is a function of radius only. In these 
lenses, the index of refraction is varied so that the radiation from a 
point source placed at the edge of the lens produces a collimated beam 
on the opposite side of the lens. The index of refraction is 1.0 at the 
edge of the lens so that it is matched to free space. 

Spherical or cylindrical Luneberg Lenses are constructed by bonding 
concentric dielectric shells or hoops together, or by a special mixing 
and blending process which produces a continuously varying index of re- 
fraction as a function of the radius [Ref. 2]. Cylindrical lenses are 
also produced by taking a cylinder of dielectric material and machining 
one surface until the required average index of refraction, when placed 
between two parallel plates, is achieved. 

Cylindrical Luneberg Lenses operate in either the TEM or the TE^q 
modes. If the thickness of the lens is greater than one-half wavelength 
but less than one wavelength, it operates in the TE^q mode. However, if 
the lens is several wavelengths thick it will operate in the TEM mode. 

Many radar systems and most electronic warfare systems require an- 
tennas that can scan a narrow beam over wide angles in space. As long 
as the scan rates are low, the present rotating systems are adequate. 
However, when the required scan rates are increased or the application 
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requires instantaneous repositioning of the beam, rotating antennas can- 
not fulfill the requirements due to their inertia. It can be shown that 

for some passive electronic warfare systems the optimum scan rates are in 
\ 

excess of 3000 rpm. In addition to this, the beamwidth of present ro- 
tating antenna systems are too wide to permit the desired accuracy in the 
determination of the direction of arrival of an intercepted signal. 

A new antenna system is therefore required to more fully utilize the 
capabilities of current equipment. The Luneberg Lens has the potential 
to overcome the above deficiencies and, in addition, provide a higher 
gain and reduce the weight of heavy antenna systems placed on the super- 
structure of ships. In view of this potential, this thesis was designed 
to analyze the parallel plate, continuously varying dielectric, cylin- 
drical Luneberg Lens. The analysis was to produce theoretical expressions 
for the far field radiation patterns of the lens for any feed system. 

This would provide a method of determining the performance of any size 
lens with any feed system to be investigated without going to the expense 
of buying or constructing a lens and measuring the far field radiation 
patterns. The only requirement is that the radiation patterns of the 
feed system be known. The feed system patterns may be known functions 
or they may be measured. The far field expressions were programmed on 
the IBM 360 computer to compute and plot the far field radiation patterns. 
The accuracy of this analysis was verified experimentally. 

Although the analysis of the Luneberg Lens was carried out for arbi- 
trary feed polarization, lens size, and frequency, some limits had to be 
placed upon the experimentation. The experimental limitations were: 

(1) feed systems polarized in a plane perpendicular to the plane of the 
lens, (2) the far field radiation patterns of the feeds were used in 
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calculating the far field expressions of the Luneberg Lenses, (3) all the 
experimentation was confined to x-band due to availability of equipment, 
and (4) a "roof-top” range was used for some measurements. 
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II. THEORETICAL ANALYSIS 



Electromagnetic energy propagating through a lens follows optical 
ray paths [Ref. 3], In the case of the radially-symmetric, cylindrical 
Luneberg Lens, the optical ray paths are such that the energy from an 
incident plane wave is focused into a cylindrical wave front on the sur- 
face of the lens diametrically opposite the direction of incidence as 




Figure 1. Focusing Properties of the Luneberg Lens, 
indicated in Fig. 1. 

The refraction necessary to achieve this focusing effect is obtained 
by properly varying the index of refraction. In the cylindrical Luneberg 
Lens the index of refraction, n(r) is given by 

n(r) = ll -(r/r o ) r , (1) 
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where r is any radius and r is the radius of the lens. Thus the index 

o 

of refraction varies from Si at the center of the lens to 1.0 at the 
edge. 

The theoretical analysis of the Luneberg Lens is based upon the 
principles of geometric optics from which the Ray Paths and the Aperature 
Field of the Lens are calculated and plotted. The far field radiation 
patterns are then calculated from the Aperature Field. 

A. RAY PATHS 

The generalized Snell’s Law for circular symmetry, such as is present 
in the cylindrical Luneberg Lens, where the index of refraction is a func- 
tion of the radius only [Ref. 3] is: 

n(R) sin 0 = , (2) 

where n(R) = /2-R 2 = index of refraction, 

R = normalized radius, 

0 = angle between the ray path and the outward meridian, 

C^= constant. 

The geometry of Fig % 2 and equation (2) were utilized to derive an 
integral formula'*' for the relationship between R, ¥ , and C^. 

f 1 4 2 2 

^ = C 1 J ~ [-R + 2R - C J] dR , (3) 

where ¥ is the angle between the direction of arrival of the plane wave 
and the radius through a point on the ray. A unique ray exists for each 
value of (0.0 <_ £ 1.0). 

For complete derivation see Appendix A. 
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Figure 2. Geometry for Luneberg Lens Ray Path Derivation. 



12 



Equation (3) was separated into two intervals as given by equations 

(41) and (42) of Appendix A and programmed for the IBM 360 computer. The 

ray path plots are given in Fig. 3 for ten values of C_ . Note that 3 

1 o 

was equal to and that the closest point of approach to the origin by 
a ray (r^) occured for 0 equal to 90°. Each ray is symmetric about the 
point (r = r , 0 = 90°). 

B. APERATURE FIELD 

The relationship between the Feed Distribution and the Aperature 
Field of the Luneberg Lens was derived utilizing the principles of geo- 
metric optics, for the geometry and variables defined in Figure 4. The 
energy flow into a Differential Element dgdx must be identical to the 
energy flow out of the Differential Element dxdy , since from the prin- 
ciples of Geometric Optics it has been shown in section 1I-A that energy 
flow is along Ray Paths. If a tube of rays is considered, then the energy 
flow is constant through various cross-sections of the tube. Therefore, 
the energy flow into the wedge-shaped tube, dBdx , at the feed is 
S^(B,x)d£dx , where S^(B,x) is the energy distribution per unit angle in 
the primary feed pattern. 

Equating the energy that enters the flux tube, dB dx , to the energy 

leaving the flux tube, dxdy , S (x,y) dxdy = S (x,B) dBdx , where S (x,y) 

J a 1 a 

is the energy in the lens aperature. However, S (x,y) is proportional 

a 

2 

to | E ’ | , where |E’|is the magnitude of the electric field in the aper- 

a a 

ature. Therefore, 

S a (x,y)dxdy = C 2 |E^| 2 dxdy = S^x^dgdx (4) 

where is a constant. From the geometry of Figure 4 

y = r sin 4' (5) 

o 
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Figure 3. Ray Paths Through Luneberg Lens. 
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Figure 4. Geometry for the Derivation of the Aperature Fields. 
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where r is the lens radius. r is defined as the point on a ray nearest 
o m 

the origin, and due to the symmetry of the Lens, the ray path is symmetric 
about the point r - r^ and hence ¥ = $ [Ref. 4]. Then from equation (5) 
and this fact 



dy dy ... 

d? ■ dii " r o C0s ’■ 



Substituting this into equation (4) yields 



S (x, T) 



or 



,, 1 2 = S l (x>g) dj3 = S l (? 

'a* dy C^r^cos 'J' 



| E ’ | = C /S (x, 4*) sec 

a o _L 



( 6 ) 



The Aperature field can therefore be calculated for any feed distri- 
bution for either the TEM or the TE^ mode Luneberg Lens. Two assumptions 
were made. First, that the feed energy distribution is a separable func- 
tion of x and H / . Therefore 



S 1 (x,^) = S 2 (x) S 3 (T). (7) 

The second assumption is that S^(x) and S^Cy) can be calculated or ap- 
proximated from the geometry or by computer curve fitting of a measured 
feed pattern. For either type Luneberg Lens the ( 4 /) can be equated 
to the feed distribution in the plane of the lens [Ref. 4]. S^Cx), how- 
ever, is different for each type of lens and is a function of the lens 
geometry. In the TE^ mode lens as is suggested from wave guide theory 

[S„ (x)] oc COS (^) , (8) 

z it -10 a 

where a is the thickness of the lens. In the TEM mode lens the assumption 
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was made that the S 2 ^ x ^ constant over the x-dimension of the lens, and 

[S 2 (x)] tem « c 4 . (9) 

The Laneberg Lens was designed to produce a constant phase in the 
plane of the lens along a line tangent to the aperature center; therefore, 
the aperature fields have a phasing 



kr (1 - cos T) . 
o 



( 10 ) 



Then for the Luneberg Lens 

E* = lE’l e j kr ° (1 - cos ^ . 
a 1 a 1 



(ID 



Combining equations (6), (7), and (9) for the TEIi mode Luneberg Lens, 



E' « /S o (40 sec T . (12) 

a j 

However, S (¥)« (E* (H-) ) 2 , so that 
j a 

E' (x,H0 = C E (40 v^ecT , (13) 

a o a 

where C is a constant of proportionality and E (T) is the magnitude of 
d a 

the feed distribution. Rewriting with the phasing information of equation 

( 10 ), 

E ' (x, 40 = C E (40 /sec 4» e^ kr ° < ' 1 _ COS (14) 

a o a 

Combining equations (6), (7), and (8) for the TE^ mode Luneberg Lens 

the aperature field is 

E' « 1 E f 1 /sTOO" sec Y cos (— ) 
a a J a 

or 

E'( X> T) - C. cos E m .Stl = ib » (1 ‘ COS ' l ' ) . (15) 

a ,b a a 
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C. ELECTROMAGNETIC FIELD CALCULATIONS 



The Electric Field of the Luneberg Lens was calculated utilizing the 
aperature field method and Huygens Principle as developed in Chapter 2 
of Ref. 5. The Electric Field equations are: 



and 



where 



E = - — V x F V X V X A, 

e wye 



H = 



A = 



F = 



— V X A — V X v X F 



4tt 



wye 
-jkR 



J(xjr )- 



4tt 



e ff - e 

' K(x',Y')- 



R 

-jkR 



R 



dS’ , 
dS' , 



J(x’ ,<?’ ) = n x h , 
a 

K(x’ ,r ) = E x n , 

3 . 

dS' = p dx’d’P’ , 



(16a) 

(16b) 

(17) 

(18) 
09) 
( 20 ) 
( 21 ) 
( 22 ) 



and the remaining variables are shown in Fig. 5. Only the Electric Field 
will be considered from this point since the Far-Field is essentially a 
plane wave, where E and H are related by the intrinsic impedance of Free- 
Space. 

2 

The Far-Field equations derived from equations (16-22) and the geo- 
metry of Fig. 5 for a feed polarized in the x-direction are: 



See Appendix B for derivation. 



18 




Figure 5. Luneberg Lens Geometry Used in the Derivation of Radiation Ex- 
pressions . 
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[1 - cos 0] [ cos $] 



(23) 



E 0 (M,O) 



= jkp £ 

An 



jkp 



E # (R,*,0) 



-jk|R-jJ -jkp cos f' 

[- [fl ][E^x\r)cos V] dx' d¥'. 



Iko e jk ° (f [ e ~j k l R ll e~ Jk0 cos 1 " ] 



4tt 



i R„ 



(24) 



[E (x 1 , Y ) ] [sin 0 sin - cos - V F I sin 0(1 + cos 0) ]dx' d^ 1 , 

3X 



and E r (R,$,0) = 0; (25) 

where E (x , ,H' ) is the magnitude of the Lens aperature field polarized 

3.X 

in the x- direct ion and 



= [R 2 + (x 1 ) 2 + p 2 + 2R(x*cos 0 cos 0 - P sin T* 



sin 0 sin 0 + p cos cos 0)]^\ 



(26) 



Tlie far-field equations for a feed polarized in the y-direction are: 



E (R,0,Q) = 



jkp e 



jkp 



4tt 



-jk|R]J -jkP cos 

[S rh 1 



(27) 



[E (x',T')][cos sin 0 (1+cos© ) - sin 0 sin '{''Jdx' dY' 
ay 



E (R,4>,0) = 



_ jkP e 



jkP 



An 



[cos $ (cos 0 + 1) ] 



(28) 



[- 



rJ k l R il 0 ~3 k P cos T ' 



-][E (x'^'Jcos T' ] dx' df', 
ay 



and E R (R,$,0) = 0; 



(29) 



where E (x * , Y ' ) is the magnitude of aperature field polarized in the 
ay 



y-direction and | R [ is as given in equation (26) 
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The Far-Field radiation patterns are identical as can be seen by 
comparing equations (23) and (24) to equations (28) and (27) respectively. 
Due to this identity, the Far-Field radiation patterns are independent of 
the Feed Polarization. Therefore, only one set of equations need be 
utilized in calculating and plotting the far-field radiation patterns 
for any feed polarization in the x-y plane. 

Approximating — by and | R | by |r| - |R'*a | (far field con- 

3 \\\ |r| 

dition) then equations (13) and (14) become: 

E = ^-^[1 - cos 0 ] [cos $]e^k^ P [ e CtX E (x) dx f e^ E ^cos ¥E (¥) dy 

0 4ttR J a i J a^ 

and (30) 

«. - & = jk(p - R) j * cE \« - 

[ sin 0 sin ¥ - cos 0 sin $ cos V - sin $ cos f]d^ , (31) 

where CEx - jkx sin 0 cos $ (32) 

and CER = jkP(cos V cos 0 - sin V sin 0 sin <£> - cos ¥) . (33) 

Equations (30) and (31) were programmed for the IBM-360 computer for 
calculating and plotting the normalized far-field radiation patterns. 

The computer programs are given in Appendix C, and the theoretical far- 
field radiation patterns are given in Section III. 



e CER cos YE ( Y) 



Calculations show that for the worst condition (V - 90°) at a distance 
of 2L2/X the phasing error is less than . 06A 5 (21° ) # 
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III. EXPERIMENTAL VERIFICATION 



A. DISCUSSION 

Experimental verification was undertaken to determine the validity 

of the equations that describe the far-field radiation patterns of the 

Parallel-Plate Luneberg Lenses. In order to accomplish this, two Lune- 

berg Lenses were acquired. One was a TEM mode lens, 44 inches (37.2 \ 

at 10 GHz) in diameter and 2.38 inches thick. The other was a TE^q 

mode lens, 10.25 inches (8.67X at 10 GHz) in diameter and 0.5 inches 

thick. All measurements were made with the lenses operating in the 

2 

receive mode and for far-field conditions (distance > 2L /X). Measure- 
ments on the mode lens were conducted in an anechoic chamber and 

those for the TEM mode lens were conducted in a ’’roof-top” range. The 
transmitting set-up is shown in Fig. 6, and was operated in the CW mode 
at 10 GHz. The only exception was for a frequency characteristics check 
of the TE iQ mode lens. The receiving set-up is given in Fig. 7. 

The procedure followed was to measure the far-field patterns of the 
feed and the Luneberg Lens and then using the feed data, calculate the 
far-field patterns and compare the two results. 

It was noted that in the case of all feeds with both lenses that when 
the feed was placed very close to the lens the first side lobes almost 
merged into the main beam. As the distance between the feed and lens 
increased the side lobes became distinct and dropped very rapidly. Then 
the particular placement of the feed was found which resulted in optimum 
(lowest) side lobe levels. When the separation was increased to a dis- 
tance greater than the optimum, the main beam spread out and engulfed 
the first side lobes. 
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Figure 6. Transmitter (X-Band) . 
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Figure 7. Receiving/Pattern Recording System. 





B. TE 10 MODE LUNEBERG LENS 

The TE^q Mode Luneberg Lens was fed by two different feed systems. 
The first was an "open-end” waveguide feed (hereafter referred to as 
the waveguide feed) and the other was a two element Yagi array (driven 
element plus director) resonant at 10 GHz. The far field patterns of 
the lens aperature fields for these feeds are given in Figures 8 and 16 
respectively. Fig. 9 is a comparison of the theoretical vs. the experi- 
mental far field patterns in the planes of the lens (E-plane), when fed 
by the waveguide feed at a frequency of 10 GHz. The gain over isotropic 
of this system is approximately 20 db, as compared to a 16 db standard 
horn, with a one-half power beam width of 6.5 degrees. The salient 
points of the two curves are given in Table I, below. 

TABLE I 





3 db 

Beam Width 


First Null 
Location 


First Sidelobe 
Location 


First Sidelobe 
Amplitude (db) 






left 


right 


left 


right 


left 


right 


Theoretical 


6.8° 


-8.2° 


8.3° 


-10.8° 


10.6° 


-19.6 


-19.6 


Experimental 


6.5° 


-7.4° 


O 

00 


- 9.9° 


f— 1 

o 

u> 

o 


-17.9 


-18.2 


Difference 


0.3° 


0.6° 


0.5° 


0.9° 


0.3° 


1.7 


1.4 



Figure 10 is a comparison of the theoretical and experimental H-plane 
(perpendicular to plane of lens at beam center) far field patterns for 
the waveguide feed. 

In order to determine the frequency characteristics of the TE^q mode 
lens when fed by the waveguide > feed far field patterns for 7,8,9 and 11 
GHz were also calculated and are presented in Figures 11 and 15. The 
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experimental far field patterns were different from the theoretical pat- 
terns by about the same amount as the patterns at 10 GHz. 

The experimental vs. theoretical E-plane far field patterns of the 
TEfo Mode Lens fed by the Yagi feed are given in Fig, 17. The differ- 
ences in side lobe and null locations and side lobe levels are negligible. 
The gain over isotropic of this arrangement is about 21.5 db as compared 
to the 16 db standard gain horn, with a one-half power beam width of 
6.2 degrees. 

A comparison of the lens aperature fields ( Figures 8 and 16) show 
that the Yagi array feed utilizes the lens aperature more efficiently 
than the waveguide feed. 

C. TEM MODE LUNEBERG LENS 

The theoretical and experimental E-plane far field patterns of the 
TEM Mode Luneberg Lens when fed with a dielectrically loaded sectoral 
horn feed (feed pattern and lens aperature field given in Fig. 19) are 
given in Figures 20 and 21, respectively. The difference between these 
two patterns is about 1.0 degrees in the one-half power beam width, 0.5 
degrees in the first null placement, 2.5 degrees in the first sidelobe 
placement, and 8.0 db in first sidelobe amplitudes. The gain over iso- 
tropic of this arrangement was about 25 db . The theoretical H-plane 
far field pattern is given in Fig. 22. The 3 db theoretical H-plane 
beam width is 22 degrees as compared to 25.5 degrees obtained experi- 
mentally by the Naval Electronic Laboratory Center [Ref. 2], 

When the TEM mode lens was fed by the waveguide feed (Fig. 8) the 
theoretical (Fig. 18) and the experimental far field pattern differences 
were somewhat less that they were for the sectoral horn, but due to 
range difficulties reliable experimental patterns were not obtained. 
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It is felt that the differences between the theoretical and experi- 
mental patterns were too great to warrant further investigation without 
refinement of the feed systems analysis. 
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Figure 8. Waveguide Feed Pattern and Luneberg Lens Aperature Field, 
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Figure 9. Luneberg Lens Far Field Radiation Patterns with Waveguide Feed. 
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Figure 10. Luneberg Lens (TE^ mode) H~plane Far Field Patterns for 
Waveguide Feed. 
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Figure 11. Luneberg Lens E-Plane Patterns for Waveguide Feed at 7.0 GHz. 
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Figure 12. Luneberg Lens E-Plane Patterns for Waveguide Feed at 8.0 GHz. 
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Figure 13. Luneberg Lens E-Plane Patterns for Waveguide Feed at 9.0 GHz. 
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Figure 14. Luneberg Lens E-Plane Patterns for Waveguide Feed at 10 GHz. 
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Figure 15. Luneberg Lens E-Plane Patterns for Waveguide Feed at 11 GHz . 
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Figure 16. Yagi Array Feed Pattern and Luneberg Lens Aperature Field for 
Yagi Feed. 
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Figure 17. Luneberg Lens (TE mode) E-Plane Patterns For Yagi Feed at 
10 GHz. 1U 
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Figure 18. TEM Mode Luneberg Lens Theoretical E-Plane Radiation Patterns 
With Waveguide Feed at 10 GHz. 
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Figure 19. Dielectric Loaded Sectoral Horn Feed Patterns and Aperature 
Field for Sectoral Horn Feed. 
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Figure 20. TEM Mode Luneberg Lens Theoretical E-Plane Radiation Patterns 
For Sectoral Horn Feed at 10 GHz. 
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Figure 21. TEM Mode Luneberg Lens Experimental E-Plane Pattern for 
Sectoral Horn Feed at 10 GHz. 
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Figure 22. TEM Mode Luneberg Lens H-Plane Radiation Patterns for Sectoral 
Horn Feed at 10 GHz. 
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IV. CONCLUSIONS AND RECOMMENDATIONS 



The results of the experimentation did not show conclusively the 
approach of this thesis to the problem of predicting the far field radia- 
tion patterns of the Luneberg Lens is valid. However, the results did 
indicate that the method is valid for the TE^q mode lens. There are at 
least two reasons for the differences, assuming that the approach is 
valid. 

First, the far field patterns of the feed systems were used in cal- 
culating the aperature fields and far-field patterns of the lens, when 
in fact, the lens was not in the far field of the feeds except for the 
case of the Yagi array on the TE ^ lens. This resulted in errors in 
the aperature fields used in calculating the far field patterns. 

The results of the TE^q mode lens seems to confirm that this was in 
fact one of the major sources of error. The far field of the waveguide 
feed is approximately 4 cm. and it was operated at 1. 7 cm from the lens 
for optimum side lobe levels. The differences between the theoretical 
and experimental patterns were an average of 0.6 degrees in first side- 
lobe placement and 1.5 db in first sidelobe amplitude. However, when the 
Yagi array feed, with a far field distance of 0 . 4 cm., was used at a dis- 
tance of 1.7 cm from the lens, the differences in first sidelobe place- 
ment and amplitude were negligible. 

The far field of the dielectric loaded sectoral horn is approximately 
9.6 cm. It was operated at about 2 cm from the TEM mode lens and the 
results were different by about 8 db in the first sidelobe levels and 
2.0 degrees in the first sidelobe placement. When the TEM mode lens 
was fed with the waveguide feed the first sidelobe differences were on 
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the order of 4.0 db in amplitude and 1.0 degrees in placement. The far 
field patterns of the TEM mode lens with the waveguide feed were not pre- 
sented in section III due to the lack of confidence in the results. 

Second , the assumption was made that the energy distribution over the 
x-dimension of the TEM mode lens is constant, while in fact it cannot be 
due to the metal plates at the top and bottom of the lens. An attempt 
to measure the aperature field distribution over the x-dimension of the 
lens failed due to the disturbance of the field caused by the dipole 
probe used. The x-variation of the aperature field did, however, appear 
to be a standing wave between the two metal plates. 

It is the opinion of the author that the approach of this thesis to 
the far field analysis of the Luneberg Lens is valid, and that the dif- 
ferences between the theoretical and experimental results is almost com- 
pletely due to the feed pattern information used in calculating the 
theoretical patterns. 

In line with the above thoughts the following recommendations for 
future investigations are made: 

A. That the near field (theoretical or experimental) of the feeds 
be used in the calculation of the far field radiation patterns of the 
Luneberg Lenses. 

B. That a broad-band Yagi array, such as Yagi constructed from 
capacitative lumped loaded cylindrical elements, be used to feed the 
Luneberg Lens. 

C. . That multiple Yagi arrays be used to determine the possibility 
of using electronically scanned Yagi array (such as in item B) to feed 
the Luneberg Lens in a direction finding system. 
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D. That the computer programs of Appendix C be modified to calculate 
the far field patterns of the Constant-K Lens. This would require a 
change of function subprograms CFCT 3 and CFCT 4 to reflect the geomet- 
ric optics of the Constant-K Lens. 
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APPENDIX A 



RAY PATH CALCULATIONS 
A. DERIVATION OF RAY PATH FORMULA 

The generalized Snell’s Law for circular symmetry is given in Ref. 3 
as 

n(R)R sin 0 = C (34) 

when 

n(R) = Index of refraction 
R = Normalized radius 

0 = Angle between ray path and outward meridian 

C^ = Constant 

From the geometry of Figure 23, 



sin 0 = 



RdT 



R 



/(RdY) 2 + (dR) 2 /r 2 + S 



cm. 2 

dT ; 



(35) 



Substituting equation (35) into (34) yields 



c i = 



nR 



L 2 , dR. 2 
/R + ( df } 



or dH* = 



C 1 dR 



r/l 2 R 2 - C 2 



(36) 



For the Luneberg Lens (normalized radius), 



n = A - ~ 2 



R 



(37) 



which yields 



dV = 



(^dR 



r/-r 4 + 2R 2 - C 2 



(38) 



or 



T = C., 



h-R 4 + 2R 2 - C 2 ] 1/2 dR. 
K I 



(39) 
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Figure 23. Geometry For Luneberg Lens Ray Path Derivations. 
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Integrating this equation 



V = j sin 1 A-C* - j sin 1 [(R 2 -C 1 2 )/(R 2 /l-C 1 2 )], 



(40) 



This equation gives four values of ¥ for each value of and R and 
hence two rays. Geometrical considerations produce the following for- 
mulas : 

y n = - sin” 1 /l-C 2 + Y • (41) 



'1 2 



and 



Y 2 = it — Y , 



(42) 



where 



R 

T 'I 



J r, HZ - 2 : 2 

1 R /-F 



dR , 



-R + 2R - Ch 



y,, corresponds to V for y < ¥ < 77 , 

1 m — — 



H* corresponds to ¥ for T < V < T 

2 o m, 



and 



y is the value of T for which R = r . 
m m 



B. CALCULATIONS OF r . 

m 

For each ray path the closest point of approach to the origin, r , 
occurs for 0 = 90°. Therefore, from equations (34) and (37) 



r [n(r ) ] sin 90° = C . 
mm 1* 



or 



r v4-r 2 = C . 

mm 1 



Squaring both sides yields 



4 2 2 

r - 2r + C- = 0 

m m 1 
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or, by the quadratic formula, 



r = 1 ± /l-C 2 . 
m 1 

The negative square root must be taken since the positive square root 
gives an r^ greater than the normalized radius. Therefore, 

r = (1 - /l-C 2 ) 1/2 , ' (43) 

m 1 

which is the closest point of approach of the ray path to the center of 
the lens . 
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APPENDIX B 



ELECTROMAGNETIC FIELD DERIVATION 

A. DERIVATION OF ELECTRIC FIELD EQUATIONS. 

The field equations for any aperature on which the electric and mag- 
netic currents are known can be calculated utilizing the Aperature Field 
method and Huygen's Principle as developed in chapter 2 of Ref. 5. The 
electromagnetic field equations for any aperature [Ref. 5], and the geo- 
metry of Figure 24 are: 



and 



where 



E = - — V x F — V x V x A 

e mye 



H = 



A « 



— V X H v X V X F; 

p wye 



— f 

4tt J _ 



j(x',r)e jk l R il 



dS ’ 



R, 



(44) 

(45) 

(46) 



F = 




K(x' > 4 > 1 )e 



-jk|Ri| 




dS' 



J(x' ,¥') = n x H 
. a 



K(x' ,4'' ) = E x n 



k 



2tt 

X 



(47) 

(48) 

(49) 

(50) 



dS' = P dx' dY' (51) 

2 2 1/2 
|R^| = [R +(R' ) -2R(x , cos < ^ > cos0-psin4' , sinQsin , M-pcos4' , cosO) ] , 



n = -sin -'i' * a + cos 4 r, a = unit normal to aperature. 

y z 



(52) 

(53) 
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Figure 24. Luneberg Lens Geometry Used in the Derivation of Radiation 
Expressions . 
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Equation (1) was divided into two parts: 



where 



and 



E = E 1 + E 2 



E = - - V x F 
1 e 



E. — V x V x A, 

2 wye 



and considered separately. 

From equations (47), (49) and (54) 



(53) 

(54) 

(55) 



E = - - V x [-£- 

1 e 4 tt 



tfe'-Jj.) e ~3 k l R ll j S . ). 



(56) 



R„ 



Taking all constants outside the curl equations and interchanging the 
order of the integration and curl operators, 



E i --47 |{ v * I K(X,, | rJ Jk|RJ l ds ' 



(57) 



Using the vector identity 



V X Kf = Vf X K + f(V X K) 



(58) 



where 




equation (57) reduces to 



E. = 



__1 

4 tt 



! e-JklRql _ _ 

(jk + TTTt) (-- r D -T-)(a r x K(x',y')) dS' 



1^1 ~Tr]1 



since 



Vf - (^- — i_> e -j fc l R ll 
I Rll I Rjl 



and 



V x K = 0. 



(59) 



( 60 ) 



K is independent of the point of observation. 
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From equations (46) and (55) 



E — V x V x [-7^- 

^ wpe 4 tt 



J (x* ' )e jk l ^ 



ds’]. 



R„ 



Taking all constants outside the curl equations and interchanging the 
order of integration and the curl operators 



E = — 2 — 
2 4Tro3e 



— — T(x T 1 )e 

V X v X ] ds*. 



(61) 



Applying the vector identity of equation (58) to the integrand twice 
yields 

2 

IR. I IR, I IR, I" IR, 



_ _ J(x', ¥ ')e Jk|R l' _ 

V x V x v v - \~ 



k~ . 2jk , 2 v -jk Ri r - “m 

+ — - ^ + ^> e 11 [a r x(a r ><J)] 



'll r l 

Then from equations (61) and (62) 



(62) 



E„ = 



j 



2 4 



irme 



(k 2 - ^ 



l R ll l R i 



2 e "j k l R ll _ _ 

)[a r x (a r x J)] ds’. (63) 



R, 



Equations (51), (53), (60) and (63) yield the Electric Field or any 

cylindrical aperature field at any distance from the aperature. For the 

2 3 

far field only, the 1/ |R^ | and 1/ |R^ | terms can be neglected leaving 



E _ jkP ff 
E 1 * 4, 



-j k |Ri | 

- — rg — | — (a x K(x’ , y ))dx’ d H*’ . 

I R 1 I r 



and 



m 2 t r "jk |R- I 

jk p f f e J ' 1 1 r — 



E 9 = / 

^ 4 TTCOG 



l R i 



[a^x (a^_ x J (x* , 4*’ ) ] dx’ d l/ 1 . 



(64) 



(65) 



The Luneberg Lens was designed to produce constant phase along a line 
perpendicular to the aperature center and therefore the aperature field 
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has a phase of kp(l - cosy’) so that: 



and 



E = E e jkp(1 ‘ COS V,) 
a a^ 



H = H e Jkp(1 “ COS H " ) 
a a n 



E and H are the E and H fields of the aperature. 
a i a i 

K and J can then be rewritten as 



K - (E x F)e jkP(1 - COS T,) - K,ei kP(1 - COS (66) 

a, 1 



and 



J = (n xiT )e jkP(1 " COS = J,e jkP(1 ‘ COS Y,) 



(67) 



Equations (64) and (65) can be rewritten as 



and 



where 



and 



F = J M 
1 4ir 



7 -jk| Ri| jk p (l - cos V’) _ 

F dx’ dY’ 



R ll 



ik“P II e “J k I R ll e ^ kP(1 " cos Y '> - 

f = 2 — £ 1/ £ ) F dx' d«F' 



E. = . 

2 4fo)£ 



R, 



F. = a x K(x',il>') = a x (E x n) , 
1 r ’ r a^ ’ 



F„ = a x (a x J (x'.t') = a x [a x (n x H )] 
2 r r 1 * r r a^ 



( 68 ) 



(69) 



(70) 



(71) 



From equation (70) and vector identity 



A x (B x C) = (A • C)B - (A • B)C, 



(72) 



or 



where 



F = (a ♦ n)E - (a • E )n, 
1 r a^ r a^ 



F 1 C 7 E ai C 8 n » 



( 73 ) 



C = a • n, 
7 r 
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and 



C q a * E 
8 r a, 



( 75 ) 



Assuming that the aperature field is polarized in the x-direction (i.e. 



E - E a and H = H a ) 



a^ ax x 



a ± ay y 



C-j = - sin 0sin $sin + cos 0cos V' 



(76) 



and 



C 0 = E sin 0cos 
8 ax 



(77) 



Combining equations (73), (76) and (77) gives 



F- = (-sin 0 sin $ sin + cos 0 cos ¥)E a 
1 ax x 



(78) 



+ E sin 0 cos $ sin V' a 
ax y 

- E sin 0 cos $ cos a 
ax z 



Converting to spherical coordinates 



where 



F- - F a + F a + F a 
1 lr r 10 0 $ 



F lr=° . 



F,_ = E cos cos $ , 
10 ax 



(79) 

(80) 
(81) 



and 



F, , = E (sin 0 sin W’ - cos 0 cos ¥ ’ sin $) 
1$ ax 



(82) 



When the above equations are substituted into equation (68) 



E = 

Air 



| R 1 1 jlcp (1 - cos V) 



I R„ 



[F- a + F a + F, a )dx f dY 1 (83) 
lr r 10 0 1$ $ 



or 



where 



E- + E a + E a + E a 
1- lr r 10 0 1$ $ 



E. = 0, 
lr 



(84) 

(85) 
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E = 

1$ 4tt 



e ~jk| RjJ g jkP (1 - cos y') 



R, 



F_ dx’ dT 
1 $ 



( 86 ) 



and 



j = isp 

10 4tt 



^-jk| R]J e jkp (1 - cos y f ) 



F 10 dx ' d ' t ' 



(87) 



Similarly, from equations (71) and (72) 



F 2 - a r , (C 9 n - C 7 H a ) 



or 



where 



F 2 = C 9 U r X n) ' C 7 (a r X V > 



C Q E a • H 
9 r a 



( 88 ) 



(89) 



and 



a x n = (sin 0 sin $ cos y % + cos 0 sin y')a 
r x 

- sin 0 cos $ cos Y 1 a 



- sin 0 cos $ sin a 

z 



(90) 



Again assuming an aperature field polarized in the x-direction 



H = H a , 
a 1 ay y 



a x H = - H cos 0a + H sin 0 cos $ a , 
r a ay x ay z 



(91) 

(92) 



Cl = H sin 0 sin $ , 
9 ay 



(93) 



and C is as given by equation (76). From equations (76) and (88-93) 
7 / 

— 2 2 2 — 

F n = H (sin 0 sin $ cos Y* + cos 0 cos V')a 
2 ay x 

. 2 

- H sin 0 sin $ cos $ cos Y a 
ay y 



- H sin 0 cos $ cos 0 cos a , 
ay z 



(94) 



Converting to spherical coordinates 



F 2 F 2r a r + F 20 3 0 + F 2* a * 



(93) 
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or 



7 = J k p 
2 4ttco£ 



where 



ii 



,-jk| Rj £ jkp (1 - cos ¥') 



R, 



< F 2r a r + F 20 a 0 + F 2»V dx ' dT ' <96) 



and 



F 2r - °> 



F. = H cos 0 cos Y'cos $ } 
20 ay 



F„ = -H sin $ cos ¥ ' 
2$ ay 



(97) 

(98) 

(99) 



When the above equations are substituted into equation (69) 

2 n tf _ e - j k l R lle jkp(1 " COS 41 "> 



7 = 3 k P 

2 47fcoe 



or 



M/ 



i- 



)[F 2r a r + F 20 a 0 + F 2t a » ldx ’ d4 " 



where 



and 



E 2 = E 2r a r + E 20 a 0 + E 2$V 



E 2r = 

•v 2 

E = J - k p 
20 4iuoe 



7 = J k ..P 

J 2$ 4Trtoe 



,-j k | kj e jkp (1 - cos ¥') 

\\\ 

, - j k | Rj_l J k P (1 ~ cos 



F 20 dx’df 



F oa . dx ’ dy ' 
2 $ 



( 100 ) 

( 101 ) 

( 102 ) 

(103) 

(104) 



From equation (53) 



where 



and 



E = E + E = E a + E a + E a A 

1 2 rr 00 $ <I> 



E = E + E. 
r lr 2r 



E 0 = E 10 + E 20- 



E = E 1 + E . 
$ 2$ 



(105) 

(106) 

(107) 

(108) 



Substituting equations (85) and (102) into (106), (87) and (103) into 
(107), and (86) and (104) into (108), 

E = 0, 
r 
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E = ( 
0 4tt J 



-jk| R l| jkp(i - COS y ' ) 



[- 



R, 



-1 [E cos y 1 cos $ + 
ax 



and 



— H cos 0 cos y'cos $ ] dx 1 dy ! , 
we ay 



f = 

$ 4tt 



fj 



t- 



r jk|Ri| jkp(i - cos y') 



] [E (sin 0 sin y ' - 
ax 



- cos 0 cos y’sin $) H H sin $ cos y ! ] dx 1 dy ! , 

ojc ay 

Considering the far field relationship of E and H at the aperature of 
the lens 

H =/^ E 
ay y ax 

which when substituted in the above equations yields 



E = 0, 
r 5 



(109) 



and 



E = [cos $ (1 + cos 0]e 

0 4tt 



(J 



— jk | R x | jkp cos y' 



R. 



E cos dx ! dy 1 
ax 






$ Att 



R- 



ax 



( 110 ) 



[sin 0 sin - cos H f, sin $ (1 + cos o) ] dx ! dy 1 . (Ill) 



A similar derivation for an aperature field polarized in the y-dir- 
ection yields 



E = 0, 
r 



( 112 ) 



jkp jkP If e-i k l E ].le-J k » cos 



E = “ e~ 

0 4tt 



[E ] 
ay 



and 



[cos y'sin $ (1 + cos ©) - sin 0 sin y' ] dx' dy’, (113) 



E = 4—^- [cos 0(cos 0 + 1)] 
$ 4 tt 



G -jk I R i | e _ 3 k d cos y* 



[E cos y'] dx' dy' , 



ay 



(114) 
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where is the magnitude of the aperature field polarized in the y- 

direction. 

The two sets of equations (109-111) and (112-114) when combined will 
produce the Far-Field expressions for any aperature field polarized in 
the x-y plane. If one is interested in the Far-Field patterns only, 
then comparison of equations (110) and (111) to equations (113) and (114), 
respectively, indicates that the normalized Far-Field patterns are iden- 
tical for any aperature field polarized in the x-y plane. Therefore, 
in order to obtain the Far-Field patterns, only one set of the equations 
need be considered. 
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appendix c 



LUNEBERG LENS RAY PATES 



IMPLICIT REAL*3(A-H,R,P) 

VARIABLES USED IN THE PROGRAM ARE: 

CI-RAY PATH CONSTANT. 

RMIN-CLOSEST POINT OF APPROACH OF RAV PATHS TO 
LENS CENTER. 

YR AY- Y-COMPONENT OF A RAY PATH POINT. 

XRAY- X-COMPCNENT OF A RAY PATH POINT. 

R- NORMALIZED LENS RADIUS. 

REAL LABEL /' </ 

R E A L* 8 TITLEU2) 

READ TITLE TO BE TYPED AT BOTTOM CF GRAPH 
RE AD ( 5 » 100) TITLE 
00 FORMAT ( 6A8 ) 

FCT5 IS A FUNCTION WHICH WHEN INTEGRATED DESCRIBES 
THE RAY PATHS THROUGH THE NORMALIZED LUNEBERG LENS. 
EXTERNAL FCT5 

DIMENSION XRAY ( 1 100) , YRAYI l 100 ) , AUXI25) 

COMMON CL 

INITILIZE CL FOR THE L ST RAY. 

C 1=0. 9D0 

PI =3. 14 159265 DO 

ENTER LOOP TO OBTAIN DIFFERENT RAYS 
DC L 1=1,9 
CALCULATE R M I N 

RMIN=DSQRT ( 1 . ODO-D SQRT ( L . OD 0-C 1**2) ) 

R= 1 .000 

ASTP = 0. 002 5D0 

N UM= ( ( R-RM IN ) /ASTP ) + l 

RHC=D ARS I N ( C 1 ) 

N=2*NUM 
K=N 
M = L 

PH=PI/2.0D0-DARSIN(DS0RT( l.0D0-Cl#*2) ) 

RR = R 

X 1 = R* DC CS ( RHO ) 

I F ( L.O. LE. XL) GO TO 2 

CALCULATION OF RAY EXTENDING IN FRONT OF LENS. 

K = K+L 

X L=Xl + 0 .005 
X R AY ( K ) =X 1 
YRAYI K) =C1 
GO TO 3 

CALL INTEGRATING SUBROUTINE TO SOLVE INTEGRAL EQ 
CALL DQATR ( RR , 1 .00,0.0000001 , 1 5 , FCT5 , DAL , IER7, AUX ) 
CALCULATE X AND Y VALUES TO BE USED IN PLOTTING RAYS 
PHI =PH+ DA L 
PH2 = P I - DA L 

XRAY IN ) =RR*DCGS( PHI ) 

YRAYI N ) =RR*DS INIPHI ) 

XRAY! M) =RR*DCOS I PH 2 ) 

Y RAY ( M ) =RR*DS I N I PH2 ) 

I F (N-M. LE. 1 ) GO TO 4 

N = N— 1 

M=M+1 

RR = RR-A STP 
GO TO 2 

ENTER INTO CONTROL FOR PLOT IMG RAYS 
I FI I .GT .1 ) GO TO 5 
MOD= 1 
GC TO 7 
5 MQD = 2 

7 CALL DRAW! K» XRAY, YRAY, MOD, 0, LABEL, TITLE ,0,0, 4, 4, 2 , 2,8, 

16 ,0, LAST) 

IF( MOD.NE. 1 ) GO TO 300 
M 00= 2 
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On 



GO TO 400 

300 I F ( I.GE.9) GO TO 350 

GO TO 400 
350 M CD=3 

400 DO 500 J = 1 , K 
500 YRAY( J)=— YRAY(J) 

CALL DRAW ( K, XR AY , Y R AY , MOD , 0 .LABEL, TITLE, 0,0, 4, 4, 2, 2, 8, 
18, 0, LAST ) 

CL = CL-0 .IDO 
L CONTINUE 

STOP 
END 



FUNCTION FCT5 ( RAN) 

FCT5 IS A CUMCTIUM 'WHICH WHEN INTEGRATED DESCRIBES 
THE RAY PATHS THROUGH THE LUNEBERG LENS. 

IMPLICIT REAL *8(A-H,R,P) 

COMMON Cl 

FCT5=C1/(RAN*DSQRT(-(RAN**4 ) +2 . 0 DO* ( RAN**2 ) - ( C 1=^*2 )) ) 

RETURN 

END 
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FEED PATTERN FITTING 
AND 

APERATURE FIELD CALCULATIONS 

C THIS PROGRAM WAS DESIGNED TO FIT A LEAST SCLARE 

C PCLYNOMINISL TO A NUMBER OF POINTS FOR A MEASURED 

C FEED PATTERN AND THEN PLOT THE FEED FUNCTION BY USING 

C THREE TIMES AS MANY POINTS (K2) IN ORDER TC DETECT ANY 

C EXCURSION THAT THE PCLYNOMINISL FIT MAY HAVE TAKEN 

C BETWEEN TWO A DJ AC I ENT POINTS . 

IMPLICIT REAL*8 ( D) 

DIMENSION DTHE T A ( LOO) , DM AG ( 100 ) , D A ! 20 ) , DW I ( 80 ) , DY ( 80 ) , 
IDELY! 80 ) ,A(20) ,DSA( 20) ,THETA( 200) , EMAG ( 20C ), DTHE! 100) 

2 , EMA( 20 0 ) 

REALMS DT I TL ( 10)/10*« '/ 

REAL LABEL/' '/ 

REAL*8 T I TLE ( 12 ) 

C READ GRAPH TITLE 

P E AD ( 5 i 100) TITLE 
100 FORMAT ( 6A8 ) 

C M IS THE DEGREE OF PCLYNOMINIAL TC eE FITTEC 

C N IS THE NUMBER OF POINTS IN THE CURVE TO BE FITTED. 

READ! 5,20C)M,N 
200 FORMAT ( 13/ 12) 

C READ IN DATA PO I NT S-DT H ET A IS THE ABSISSIA VALUE AND 

C DM AG IS THE SQUARE OF THE ORDINATE VALUE-SCUARE LAW 

C DETECTOR WAS USED TO OBTAIN DATA. 

READ! 5, 400) ( DTHET A ( I ) , DMAG ( I ) , 1 = 1, N) 

DC 60 J = 1 , N 

60 DM AG( J ) =DSQRT ( DM AG ( J ) ) 

400 FORMAT! F10 .7, FL2 .6 ) 

DO 50 I = 1 1 N 

50 DTHE ( I ) =DT HET A ( I)*(2. 14159 2 65D0/180.D0) 

DO 5 1=1 ,N 

5 CW I ( I ) = l .000 

C LSCPL2 WILL CALCULATE THE LEASE SQUARE POLUNOMINIAL 

C CCEFFC I ENTS . 

CALL LSQPL2(N,M, DTHE , DM AG , D W I , DY , DEL Y , DA , D SA , D T I TL ) 
L=l+I ABS(M) 

DO 6 I = 1 » L 

6 A ( I ) = DA ( I ) 

WBITE(6,300) ( I , A! I ) , 1 = 1, L ) 

300 FORMAT! • • , ' A ( • , I 4 , • ) = • , F 1 8 . 6 ) 

TFET=DTHE ( 1 ) 

K 2 = 3* N 

C CALCULATE SIZE OF STEPS TO BE USED IN PLOTTING CURVE. 

STEP=( DTHE ( N ) -DTHE ( 1) )/K2 
N R I TE ( 6 , 3 02 ) 

302 FORMAT!' • , 4X , • THE T ' , 1 2X , ' E MAG ' ) 

DC 3 K=1,K2 

th:eta(K)=thet 

EM = A ( 1 ) 

DC 4 J = 2 , L 

4 EM=EM+A( J)*THETA(K)*»(J-1) 

EMAG(K)=EM 

WP I TE ( 6 , 301 ) THE T , EM AG ( K ) 

301 FORMAT! • • , F 10. 5 , 5X ,F l G .7 ) 

TFET=THET+STEP 

3 CONTINUE 

C CALCULATE THE APERATURE FIELD STRENGTH 

DC 7 K= 1 , K2 

EMA(K) =EM AG! K) *SQRT ( L.O/COS! THETA (K) )) 
WRITE(6,301)THETA(K) ,EMA(K) 

7 CONTINUE 

DO 8 K= 1 , K 2 

THETA (K )=THETA(K) *180.0/3. 1415 

8 CONTINUE 

CALL DR AW ( K2 , TH ET A , EM AG , 1 , 0 , L A E E L , T I TL E , 2 C . 0 , 0 , 0 , 3 , 2 , 

12.6.7. 1, LAST) 

CALL DRAW ( K2 , THETA, EMA ,3 ,0, LABEL, T ITLE, 2C . C , 0 , 0 , 3, 2, 

12.6.7. 1 , LAST) 

STOP 
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END 
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PH E-COMPCNENT LUNEbERG LENS PATTERNS 



C THIS PROGRAM WAS DESIGNED TO CALCULATE P h E-COM PON EN T 

C OF THE RADIATION PATTERNS OF THE CYLINDRICAL LUNEBERG 

C LENS WITH BEAM CENTERED AT O-DEGRE ES , FOR A VERTICALLY 

C POLARIZED FEED. THIS PATTERN IS ICENTICAL FCP ANY FEED 

C POLARIZATION IN THE X-Y PLANE 

IMPLICIT COMPLEX ( C ) 

COMPLEX CJ/tO. 0,1.0)/ 

M = C 

C STEP IS THE STEP-SIZE IN DEGREES BETWEEN POINTS WHERE 

C THE FIELD STRENGTH IS CALCULATED. 

C THETAI IS THE INITIAL ANGLE IN DEGREES FOR WHICH THE 

C FIELD STRENGTH IS CALCULATED. 

C ThETAF IS THE FINAL ANGLE IN DEGREES FOR WHICH THE 

C FIELD STRENGTH IS CALCULATED. 

C ANGPD IS THE ANGLE IN CEGREES THAT CONTROLS THE PLANE 

C IN WHICH THE LENS PATTERNS ARE CALCULATED, IF 90 

C CEGREES THE THE E-PLANE PATTERNS ARE CALCULATED, IF 

C O-DEGREES THEN THE VERTICAL BEAM PATTERN IS CALCULAT- 

C EC* 

R E AD ( 5 , 2 ) STEP, THETAI, THETAF, ANGPC 
2 FORMAT (FL2.3/FL2.3/FI2.3/FI2.3) 

PI E = 3. 1415926535 
ANGP=ANGPD“( P IE/ 1 80. 0) 

C CALCULATE NUMBER OF STEPS ,NUM, WHICH IS ALSO THE 

C NUMBER OF POINTS TO EE PLOTTED BY SUBROUTINE CRAW. 

NUM=( THETAF-THETA I ) /STEP + 1 

C MG IS THE INTEGER THAT CONTROLS THE NUMBER OF PLOTS 

C TO BE PLOTTED BY DRAW IN SUBROUTINE EPhEE. 

C RADIUS IS THE RADIUS OF THE LUNEBERG LENS. 

C FREQ IS THE FREQUENCY FOR WHICH THE PATTERNS ARE 

C CALCULATED. 

C THICK IS THE HEIGHT OR X-DIMENSICN OF THE LENS. 

12 READ(5.4) MO , RADIUS, FREQ, THICK 

A FORMAT! I 1/F L2.3/F L2. 3/F 12. 3) 

I F ( MO . EQ . 0 ) GO TO 14 
N = N+ 1 

LA WRITE(6,5)THETAI, THETAF, STEP, RAD I US , FR EC , NON , MO , 
1THICK, ANGPD 

5 FORMAT ( ' ', 'THETAI*', Fll. 6,/,' ' , • Th ET A F = • , F 1 1 .6 , / , ' 

l ' , 'STEP* • ,F1 1.6, / , • ', 'RADIUS*' ,F11. 6,/,' ', 

2* FREQUENCY*' ,F11 .6,/,' ' , ' NUM= ' , 16 , / , ' ','MC ',16,/' 

3, 'THICK*' ,F1L.6,/, ' ', ' PHE = ' , F 1 I . 6 ) 

R = 200 . 0 

C AK IS THE PROPAGATION CONSTANT. 

AK=20.0*PIE*FREQ/3.0 

C JKR=(C J*AK-*RADIUS/ (4.0*P I E*R) )*CEXP(CJ*AK=* (RADIUS-R ) ) 
C CJKRN IS THE NORMALIZED COEFFICIENTS 

CJKRN=CJKR/CABS( CJKR) 

WRITE! 6 ,100)CJKRN 

100 FORMAT!' ' , • C JKRN= • , F 1 l . 9 , 2X , • J ' , F 1 1 . 9 ) 

CALL EPHEE! AK, THICK, ANGP, THETAI , STEP, NUM , C JKRN , RAC I US 
l , MC , FR EQ , M ) 

C EPHEE WILL CALCULATE AND PLOT THE PHE-CO M PC N E NT S 

C OF THE E-FIELD OF THE LUNEBURG LENS 

I F ( MO . EQ . C ) GO TO 13 

13 STOP 
ENC 



SUBROUTINE EPHEE ( AK 1 , TH I , ANGP l , ANGT I , S T E , N ON , C JK R N , 
IRA , MOD, F R E , M 1 ) 

C THIS SUBROUTINE WAS WRITTEN TO CALCUUATE AND PLOT 

C THE PHE-COMPONENT OF THE ELECTRIC FIELD OF THE LUNE- 

C BERG LENS FOR ANY FEED FUNCTION OR ANY MEASURED FEED 

C DISTRIBUTION THAT CAN BE APPROXIMATED BY A CURVE 

C FITTING TECHNIQUE. 

IMPLICIT COMPLEX(C) 

C THE NEXT FOUR STATEMENTS CONTROL THE LABELING AND 
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C TITLING OF THE GRAPHS PLOTTED BY SUBROUTINE CRAW. 

REAL LABEL /.' '/ 

REAL *3 riTLt(L2) 

READ! 5 , 100 ) TITLE 
LOO FORMAT ( 6A3 ) 

COMPLEX C J/( 0 .0, 1 .0)/ 

EXTERNAL CFCTL,CFCT3 

C CFCT1 IS A FUNCTION DESCRIBING THE FAR FIELD OF THE 

C LUNEBERG LENS IN TERMS OF THE APPERATURE FIELC 

C DISTRIBUTION IN THE X-D I RECT I ON . 

C CFCT3 IS A FUNCTION DESCRIBING TEE FAR FIELC OF THE 

C LUNEBERG LENS IN TERMS OF THE APPERATURE FIELD 

C DISTRIBUTION AS A FUNCTION OF THE 

C ANGULAR DISPLACEMENT FROM BEAM CENTER. 

DIMENSION THETA! 40 0 .CEPHE (400) , EPHEM (400 I , R AD I (5 ) , 
1FREQ(5) 

COMMON/ A L I M/ AL I ML , ALIM2.THIC 
COMMON AK, ANGP, THETR,RAD,L3 ,LL ■ 

T H I C=T H I 
R AD = R A 
RADI! ML I = R A 
FREQ( ML ) = FRE 
A K = AK1 
ANGP=ANGP l 
PIE=3.L4L59265358 
THET=ANGT I 
TH2=-TH 1/2.0 
TH3=T H I / 2 . 0 
A N C RM = 0 . 0 

C A NORM IS A VARIABLE USED TO NORMALIZE THE ELECTRIC 

C FIELD STRENGTH. 

C ENTER LOOP TO CALCULATE THE E-FIELO FOR NUM-POINTS 

DO L 1 = 1, NUM 
Li= I 
L 3 = 1 

THETR = THET *( P IE/ 130.0 I 

CALL CWEDF(CFCTL,TH2,TH3, L5,CEPH1) 

L 1 = 2 

900 X L = AL I M 1 

XU = AL I M2 

CALL CWEDF (CFCT3 , XL , XU ,375,CEPH3) 

IFIL3.NE.-U GO TO 1000 
L 3 = 2 

GO TO 900 

1000 CEPH=CEPHl*CEPH3 

IF(ANORM.GE.CABS (CEPH) ) GO TO 2 
AN0RM=C A3S (C-EPH ) 

2 THETA! I ) = T HET 
CEPHEU ) =C EPH 
THET=THET+STE 

L CONTINUE 

C ENTER L COP TO NORMALIZE AND CALCULATE THE MAGNITUDE 

C OF THE E-FIELD (EPHEM). 

DC 3 J= 1 , NUM 

CEP HE l J ) =CEPHE ( J ) /ANORM 

EPHEM (J l=CABS! CEPHE ! J I ) 

3 CONTINUE 

WRITE !6,34 I RAD I (ML) , FREQ! ML) 

84 FORMAT (• L ', • PHE- COMPONENT OF THE E-FIELD FCR LENS', 

L ' RADI US OF* , F3.4 ,2X, ' METERS AT ' , F8 . 3 , 2X , ' GHZ ' , // ' • , 6X 

2, 'THETA ' , L2X, *F I ELD STRENG TH • , I OX , • COMPL E X FIELD ', 

3* STRENGTH' ) 

J = 0 

DO 23 1=1, NUM 

WRITE (6, 85 )THET A! I ) , EPHEM! I ), CEPHE! I ) 

85 FORMAT!' • , F 1 L . 3 , 8 X , F 14 . 6 , L 2X , F 1 0 . 6 , F L 0 . 6 ) 

J=J + 1 

IFtJ.LT.72 ) GO TO 23 
WRITE ( 6,84)RAD,FRE 
J =0 

23 CONTINUE 

C ENTER LOOP TO CALCULATE THE FIELC MAGNI TUCES IN DECI- 
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C BELS. 

DG 20 I =1 , NUM 

EPHEMU ) = 2 0. 0*AL0G 1 0 ( EPHEM ( I ) ) 

20 CCNT IN'JE 

CALL OR AW (NUM, THETA, E PHEM , MOD, 0 , LAE EL, TITLE, 0,8. 0,7, 2, 
12,2,6,3, l, LAST) 

RETURN 
END 



SUBROUTINE CWEOF ( C F , XL , XU , NX , C AN S ) 

CWEDF WAS OBTAINEO FROM LCDR . E.G. NEELY . 

CWEOF IS A SUBROUTINE WHICH WILL NUMERICALLY 
INTEGRATE A USER SUPPLIED FUNCTION BETWEEN SPECIFIED 
LIMITS. (SINGLE PRECISION) 

CF - NAME OF FUNCTION SUBPROGRAM. MUST BE 
LISTED IN AN EXTERNAL STATEMENT. 

XL - LOWER LIMIT OF INTEGRATION 

XU - UPPER LIMIT OF INTEGRATION 

NX -APPROXMATE NUMBER OF NGCES AT WHICH TO 
EVALUATE THE FUNCTION. 

CANS - RESULT OF INTEGRATION 

IMPLlCI T COMPLEX--^ ( C) 

REALMS DXDX , XX 

REAL*4 CW (6 )/3 2 . ,216. ,27. ,272 ., 27., 216./ 

COMMON AK,ANGP, THETR ,RAD, L3 ,L1 
IF(NX.LE.O) GO TO 900 
N=( ( NX + 4 ) / 6 ) *6+1 
DX=( XU--XL) /FLOAT ( N- 1 ) 

DXCX=DBLE( DXJ 
NWIX=N/6 
X= XL 

CANS=-CF(X )*41 .0 
IF(L3.EC.-1) GO TO 1000 
IF(Ll.EQ.l) GO TO 500 
I F ( L3 . N E . 1 ) GO TO 500 
L 2 = 2 

DO 800 MX= 1 , NKI X 
DC 700 KX = 1 , 6 
CANS=CANS+CW(KX)*CF(X) 

XX=DBL E (X ) 

X=SNGL( XX+DXDX ) 

CONTINUE 

CONTINUE 

CANS=(CANS+41 .0*0 F ( X ) ) *DX/ 1 40 .0 
RETURN 

WRITE (6, 901) N 

FORMAT ( • OERROR IN CALLING PARAMETER ***** N = ' , 

1 I 5 , • ****** / / ) 

RETURN 

END 



FUNCTION CFCT1 (X ) 

C CFCT1 IS A FUNCTION WHICH DESCRIBES THE FAR-FIELD 

C VARIATIONS OF THE E-FIELD AS A FUNCTION OF THE LUNE- 

C BERG LENS HEIGHT. 

IMPLICIT COMPLEX (C ) 

COMPLEX CJ/IO.O, 1.0)/ 

COMMON AK.ANGP, THETR, RAD,L3 ,L1 
COMMON/ AL I M/ALIM1, ALI M2.THIC 
CEX=CJ*AK*SIN(THETR)*COS( ANGP)*X 
CFCT1=CEXP(CEX)*C0S(3. 14159* X/THIC) 

RETURN 

END 



90 



C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

c 



500 



700 

800 



900 

901 



1000 
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FUNCTION CFCT3IANG) 

C CFCT3 IS A FUNCTION WHICH DESCRIBES THE FAR FIELD 

C VARIATIONS OF THE E-FIELD AS A FUNCTION OF ANG 

C (THE ANGLE FROM BEAM CENTER). TO DO THIS LSCPL2 IS 

C UTILIZED TO FIT A CURVE TO THE MEASURED FEED FUNCTION 

C AND COMBINES THIS WITH THE OTHER ANGULAR VARIATIONS 

C OF THE FAR-FIELD EQUATIONS. 

IMPLICI T R E A L * 8 ( D ) 

IMPLICIT COMPLEX { C J 
COMPLEX CJ/tO.O, I .0)/ 

DIMENSION DTHETA ( LOO) , DMAG ( 100 ) , D A ( 20 ) , DW I ( 30 ) , DY ( 80 ) , 
LDELY(30),A(20),DSA( 20) , THE TA ( LOJ) , EMAG ( IOC) » DTHE ( LOO) 
REAL*6 DTITL( 10J/10*' '/ 

COMMON AK f ANGP , TH E TR , R AD , L 3 , L 1 
COMMON/ AL IM/ ALIM 1, AL I M2 , TH I C 
C L 3 WILL BE EQUAL TO 1 ON THE FIRST CALL OF CFCT3 

C AND WILL PROCEED THROUGH THE NEXT STEPS DOWN TO 600 

C IN ORDER TO CALCULATE THE POLYNOMINIAL COEF. THAT 

C DESCRIBE THE FEED FUNCTION. 

C IF L 3 IS NOT 1, THEN THE POLYNOMINIAL COEF. HAVE BEEN 

C CALCULATED BY A PREVIOUS PASS. 

I F( L3.NE. L) GO TO 600 
R EAD ( 5, 200 ) M,N 
200 FORMAT* 13/ 12 ) 

READ( 5, 400 ) ( DTHE TA ( I ) » DMAG ( I ) » I =1 » N) 

400 FORMAT { FLO .7, FL2 .6 ) 

DC 60 J = 1 » N 

60 DMAG( J )=D5QRT( DMAG( J) ) 

DO 50 I = 1 » N 

50 DTHEl I ) =D THE TA { I )*(3. 14159265 DO /180. DO) 

AL IM1=DTHE( L ) 

ALIM2=0THE (N) 

DO 5 1=1, N 

5 DW I ( I ) = 1 .ODO 

CALL LSCPL2( N,M, DTHE , DM AG , DW I , CY , DEL Y , D A , DSA , DT I TL ) 
L = 1 + 1 AB S( M) 

DC 6 1=1, L 
A ( I ) =D A { I ) 

6 THETA ( I ) =DTHETA ( I ) 

WR ITE ( 6 , 300 ) ( I , A( I ) , 1= L,L ) 

300 FORMAT { • • , ' A( ' , 14 , • ) =• , FI 8 .6 ) 

L3=-l 

RETURN 

600 E M = A ( 1 ) 

DO 4 J = 2 , L 

4 EM=EM+A(J )*ANG**( J-l) 

C E A 1 IS A FUNCTION WHICH DESCRIBES THE FEEC VARIATIONS 

C OVER THE PLANE OF THE LENS. 

E A 1 = E M 

CER=CJ*AK*RAD*(C0S(ANG)*(C0S(THETR)-1.0)-SIN(ANG)* 
1SIN(THETR ) *SIN(ANGP) ) 

C THE SQRT( 1 . 0/COS (ANG) ) TERM IS TFE GNLY PART OF THIS 

C EQUATION ASSOCIATED WITH THE LENS FUNCTION. 

C THE REMAINDER IS FROM THE APPLICATION OF THE APERATURE 

C FIELD METHOD TO THE SEMI-CIRCULAR APERATURE. TO MOD- 

C I FY THIS PROGRAM TO PREDICT THE CONST ANT- K LENS FAR 

C FIELD PATTERNS CHANGE THIS FUNCTION. THE CORRECT PHASE 

C TERMS ARE ALREADY INCLUDED. 

CFCT3=CEXP(CER)*( S IN< THETR)*SIN ( ANG)-COS( ANG)* 
1SIN(ANGP)*(C0S(THETR) +1 .0) )*SQRT(L .0/C0S(ANC-) )*EA1 
RETURN 
END 
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THETA-COMPONENT LUNEBERG LENS PATTERNS 



C THIS PROGRAM WAS DESIGNED TO CALCULATE TEE THETA CGM- 

C PONENT OE THE CYLINDRICAL LUNEBERG LENS RACIATION PAT- 

C TERNS WITH BEAM CENTERED AT O-D EGR EES , FO R ANY FEED 

C POLARIZED IN THE X-Y PLANE. 

IMPLICIT COMPLEX (C) 

COMPLEX CJ/(0.0, 1 .0)/ 

M = 0 

C STEP IS THE STEP-SIZE IN DEGREES BETWEEN PCINTS WHERE 

C THE FIELD STRENGTH IS CALCULATED. 

C THETA I IS THE INITIAL ANGLE IN DEGREES FOR WHICH THE 

C FIELD STRENGTH IS CALCULATED. 

C THETAF IS THE FINAL ANGLE IN DEGREES FOR WHICH THE 

C FIELD STRENGTH IS CALCULATED. 

C ANGPD IS THE ANGLE IN DEGREES T.HAT CONTROLS THE PLANE 

C IN WHICH THE LENS PATTERNS ARE CALCULATED, IF 90 

C DEGREES THE E-PLANE PATTERNS ARE CALCULA TED , I F 

C O-DEGREES THEN THE VERTICAL BEAM PATTERN IS CALCULAT- 

C ED. 

R E AD! 5 » 2 ) STEP.THETAI, THETAF, ANGPD 
2 FORMAT! FI 2.3 /F 12 .3 .F12.3/F12.3) 

PIE=3. 14L5926535897 
ANGP= ANGPD* ( P IE/ 18 0 .0 ) 

C CALCULATE NUMBER OF STEPS REQUIREC WHICH IS THE SAME 

C AS THE NUMBER Of- POINTS TO BE PLOTTED BY DRAW 

NUM=!THETAF-THETAI 1 /SI EP +L 

C MO IS THE INTEGER THAT CONTROLS THE NUMBER CF PLOTS 

C TO BE PLOTTED BY DRAW IN SUBROUTINE ETHETA 

C RADIUS IS THE RADIUS OF THE LUNEBERG LENS. 

C FRE3 IS THE FREQUENCY FOR WHICH THE PATTERNS ARE 

C CALCULATED. 

C THICK IS THE HEIGHT. OF THE LENS. 

L2 R E AD( 5, A ) MO , RADIUS, FREQ, THICK 

4 FORMAT! U/F12.3/FI2.3/FI2. 3) 

I F ( MO . EQ . 0) GO TO 14 

M= M+ 1 

14 WRITE(6,5)THETAI , THETAF, STEP, RAC IUS , FR EQ , NUM, MO, 
ITHICK, ANGPD 

5 FORMAT ( f • ,'THETAI=',F11.3,/, • 1 ,' THETAF = ', F U . 3 ,/ , 

1* ' ,'STEP=» ,Fii.6, /,« • , ' RADIUS*' ,FI1 .6,/, • ', 

2 ' FREQUENCY* ' , F L L .6 , / , ' ' , ' NUM= ' ,16,/,* • , • MC =,',16 

3,/' • , 'THICK*' ,FII .6,/* ' , • PFE= • , Fll .6) 

R=200. 0 

C AK IS THE PROPAGATION CONSTANT. 

AK =2 0.0* PIE* FREQ/3 .0 

CJKR=(CJ*AK*RADIUS/( 4. 0*PIE*R) ) *CE XP ( C J * AK* ( R AD I US-R ) ) 
CJKRN=CJi\R/CABS (CJKR) 

WRITE! 6 ,100)CJKPN 

100 FORMAT!' « , • C JKR N= • , F l L . 9 , 2X , • J • , F 1 1 . 9) 

CALL ETHETA! AK, THICK , ANGP, THETA I , ST EP , NUM , C JKRN , R A D I US 
IMG , F RE Q , M ) 

C EHTETA WILL CALCULATE AND PLOT THE THETA COMPONENTS 

C -FIELD PATTERNS OF THE LUNEBURG LENS 

I F ! MO .EQ. 0 J GO TO 13 
13 STOP 

END 



SUBROUT INE ETHET A ( AK 1 , TH I , ANGP 1, ANGTI , STE , NUM,C JKRN, 
1RA,M0D,FRE,M1 ) 

C THIS SUBROUTINE WAS WRITTEN TO CALCULATE ANC PLOT 

C THE THETA-COMPONENT OF THE ELECTRIC FIELD OF THE LUNE- 

C BERG LENS FOR ANY FEED FUNCTION OR ANY MEASURED FEED 

C DISTRIBUTION THAT CAN BE APPRIXIMATED BY A CURVE 

C FITTING TECHNIQUE. 

IMPLICIT COMPLEX(C) 

C THE NEXT FOUR STATEMENTS CONTROL THE LABELING AND 

C TITLING OF THE GRAPHS PLOTTED 6Y SUBROUTINE DRAW. 

REAL LAdEL /' • / 
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REAL * 8 T1 TLF( 12) 

READ(5, 100) TITLE 
100 FORMAT! 6A8) 

COMPLEX CJ/( 0.0, 1.0)/ 

C CFCT1 IS A FUNCTION DESCRIBING THE FAR PI ELD OF THE 

C LUNEBERG LENS IN TERMS OF THE APPERATURE FIELD 

C DISTRIBUTION IN THE X-DIRECTION. 

C CFCT4 IS A FUNCTION DESCRIBING THE FAR FIELD OF THE 

C LUNEBERG LENS IN TERMS OF THE APPERATURE FIELD 

C DISTRIBUTION AS A FUNCTION OF THF 

C ANGULAR DISPLACEMENT FROM BEAM CENTER. 

EXTERNAL CFCT1,CFCT4 

DIMENSION THETA (400) ,CETHE( 400 ) , ETHEM( 400 ) , RADI ( 5 ) , 

1 FREQ( 5 ) 

COMMON AK, ANGP,THETR,RAD,L3,L1 

COMMON/ ALI M/ALIML , ALI M2.THIC 

TH IC=TH I 

R AC=RA 

RACI (Ml )=R A 

F R E Q ( M 1 ) =FRE 

AK= AK 1 

ANC-P=ANGP1 

PI E=3. 141592653589 

THET= AN GT I 

T F2=-T H I / 2 .0 

TH3=TH I / 2. C 

AN0RM=0 .0 

C ENTER LOOP TO CALCULATE THE E- F I EL C FOR NUM-POINTSIN 

DO 1 1=1, NUM 
L 1=1 
L3 = I 

THETR=THET*( PIE/ 180.0) 

CALL CW EDF (CFCT 1, T H2 , TH3, 5 , C ET H 1 ) 

L 1 =2 

900 XL=AL I M 1 

XU = AL I M 2 

CALL CW EDF ( CFCT 4 , XL , XU ,375 ,CETH4) 

IF(L3.NE.-1) GO TO 1000 
L3 = 2 

GC TO 900 

1000 CETH=CETH1*CETH4 

WRITE (6 ,96 ) CETH1 ,CETH4,CETH 
96 FORMAT ( ' ' ,6E 12.6) 

C WR ITE ( 6 , L5 )CETH 

15 FORMAT ( ' • , • CETH=' ,2F10 .5 ) 

I F ( ANORM. GE . CABS ( C ETH ) ) GO TO 2 
ANCRM=CABS (CETH) 

2 THETA ( I ) =T HE T 
CETHE ( I ) = C ET H 
TF ET = T HET + ST E 

WRITE (6,95) T HET , ANCRM 
95 FORMAT ( « « ,2F12.5) 

1 CONTINUE 

C ENTER LOOP TO NORMALIZE ANO CALCULATE THE MAGNITUDE 

C OF THE E-FIELD (EPHEM). 

DC 3 J=l, NUM 

CETHE ( J ) =C E THE ( J) /ANORM 

E T H E M ( J ) = CABS(CETHE(J) ) 

3 CONTINUE 

WRITE ( 6 ,34) RADI (Ml) ,FREQ(M1) 

84 FORMAT (• 1 ', 'THETA-COMPONENT OF THE E-FIELD FOR LENS', 

1 ' RADI US OF' , F8 .4 ,2X, ’ METERS AT',F8.4,2X, 'GHZ',//' ' , 6X 

2, 'THETA ' , 12X, 'FI ELD STRENGTH' , 10X ,' COMPLEX FIELD ', 

3* STRENGTH* ) 

J = 0 

00 23 1=1, NUM 

WRITE. ( 6 , 8 5 ) T HET A ( I ) , ETHEM ( I ) , C ETHE ( I ) 

85 FORMAT ( ' ■ ,F 1 1 . 3 , 8X , F 14 .6 , 1 2X , F 10 .6 , F10 . 6 ) 

J = J + 1 

1 F ( J.LT .72 ) GO TO 23 
WRITE ( 6 , 84 ) R AD , F R E 

J = 0 
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23 CONTINUE 

C ENTER LOOP TO CALCULATE THE FIELC MAGNITUDES IN DECI- 

C BELS. 

DC 20 I =1 , NUM 

ETHEM( I 1=2 0. 0*ALOG10( ETHEMI I ) ) 

20 CONTINUE 

CALL DR AW ( NUM, TH ET A , ETHEM , MOD , 0 , LABEL, TITLE, 0,0, 7, 3, 2, 
12, 6, 8, L .LAST) 

90 RETURN 

END 



C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 



500 



700 

800 



900 
90 L 

1000 



C 

C 

c 



c 

c 

c 



SUBROUTINE CWEDF ( C F , XL , XU , NX , C AN S ) 

CWEDF WAS OBTAINED FROM LCDR. E.G. NEELY . 

CWEDF IS A SUBROUTINE WHICH WILL NUMERICALLY 
INTEGRATE A USER SUPPLIED FUNCTION BETWEEN SPECIFIED 
LIMITS. (SINGLE PRECISION) 

CF - NAME OF FUNCTION SUBPROGRAM. MUST BE 
LISTED IN AN EXTERNAL STATEMENT. 

XL - LOWER LIMIT OF INTEGRATION 

XU - UPPER LIMIT OF INTEGRATION 

NX -APPRQXMAT5 NUMBER OF NCDES AT WHICH TO 
EVALUATE THE FUNCTION. 

CANS - RESULT OF INTEGRATION 

IMPLICIT COMPLEXES (C) 

REA L* 3 DXDX.XX 

REAL* 4 CW( 6) /8 2. ,216. ,27. ,272. ,27. ,216./ 

COMMON AK, ANGP, T HE TR , R A D, L 3 , L I 
IF(NX.LE.O) GO TO 900 
N= ( (NX+4) / 6 ) *6+1 
CX= ( XU-XL ) /FLOAT ( N— 1 ) 

DXDX=D6 LE ( DX ) 

NWIX=N/6 
X = x L 

C ANS=-*C F ( X ) *4 1 . 0 
I F ( L3 . EQ •— 1 ) GO TO 1000 
I F ( LI . EC. L » GO TO 500 
I F ( L3.NE. L )G0 TO 5C0 
L 3 = 2 

DO 800 MX = 1 , N W I X 
DO 700 KX=1,6 
CANS=CANS+CW (KX)*CF(X) 

XX=DBLE (X) 

X=SNGL ( XX+DXDX) 

CCNT INUE 
CCNTINUE 

CANS=(CANS+41.0*CF(X) )*DX/140. 0 
RETURN 

WRITE (6,901) N 

FORMAT ( 'OERROR IN CALLING PARAMETER ***** N = ', 
118 ,' ****** // ) 

RETURN 

END 

COMPLEX FUNCTION CFCTl(X) 

CFCT1 IS A FUNCTION WHICH DESCRIBES THE FAR-F I ELD 
VARIATIONS OF THE E-FIELD AS A FUNCTION OF THE LUNE- 
BERG LENS HEIGHT. 

IMPLICIT COMPLEX(C) 

COMPLEX CJ/(0.0, 1.0)/ 

COMMON AK, ANGP, THETR, RAD, L3 ,L1 
COMMON/ALIM/ALIM L, ALI M2, THIC 
CEX=CJ*AK*SIN(THETR)*COS ( ANGP) *X 

C FCT1=C EXP (C EX) *C0S(X*3. 141 59/TH IC)*(1 . 0 + CC S ( THET R ) )* 
1CCS (ANGP) 

RETURN 

END 

COMPLEX FUNCTION CFCT4 ( ANG ) 

CFCT4 IS A FUNCTION WHICH DESCRIBES THE FAR FIELD 
VARIATIONS OF THE E-FI FLD AS A FUNCTION CF ANG 
(THE ANGLE FROM BEAM CENTER). TO DO THIS L $ L PL 2 IS 
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c 

c 

c 



c 

c 

c 

c 

c 

c 



200 

400 

60 

50 



5 



6 

300 



600 



4 

C 

C 

7 



C 

C 

3 

C 

C 

C 

C 



UTILIZED TO FIT A CURVE TO THE MEASURED F EE C FUNCTION 
AND COMBINES THIS WITH THE OTHER ANGULAR VARIATIONS 
OF THE FAR-FIELD EQUATIONS. 

IMPLICIT R E AL*8 ( D ) 

IMPLICIT COMPLEX { C ) 

COMPLEX CJ/IO.O, I .0)/ 

DIMENSION QTHETA! LOO) , DMAG { 100 ) , D A { 20 ) , DW I < 80 ) , OY ( 80 ) , 
IDELY{30),A(20),DSA(20) , THETA! LOO) » EMAG { LOC) ,JTHE! 100) 
REAL-8 CTITLI 10 ) / 1 0 * ' */ 

COMMON AK, ANGP,THETR,RAD,L3 ,L1 
COMMON/ AL IM/ALIM1, ALIM2.THIC 

L3 WILL BE EQUAL TC L ON THE FIRST CALL OF CFCT4 
AND WILL PROCEED THROUGH THE NEXT STEPS DCWN TO 600 
IN ORDER TO CALCULATE THE POLYNOMINI AL COEF. THAT 
DESCRIBE THE FEED FUNCTION. 

IF L 3 IS NOT 1, THEN THE POLYNOMINIAL COEF. HAVE BEEN 
CALCULATED BY A PREVIOUS PASS. 

I F { L3.NE. I ) GO TO 600 
READ( 5, 200 )M, N 
FORMAT { 13/ 12 ) 

READ! 5, 400 { DTHETAI I ) ,DMAG( I ) ,1 = 1,N) 

FORMAT ( FIO .7 , FI2 .6 ) 

DO 60 J =1 * N 

DMAGIJ )=DSQRT(DMAG(J) ) 

DC 50 1=1, N 

DTHEI I ) =D THETA! I ) *13. 1 42592 65D0/ 18 0.0 DO ) 

AL I M 1=DTHE ( L ) 

A L I M2 =DTHE ( N ) 

DO 5 1=1, N 
DW I ( I ) = l .ODO 

CALL LSGPL2I N,M, DTHE, DMAG, DW I, DY , C ELY , D A , C S A , DT I T L ) 
L=l-M 

DC 6 1= l, L 
A ( I ) =D A I I) 

THETA! I )=D THETA! I ) 

WRIT E(6, 300) ( I, A! I ) , 1= l,L ) 

FORMAT!* * ,*A(* ,14,* ) = * ,F18.6) 

L 3=- 1 
RETURN 
E M = A ( l ) 

DO 4 J = 2 , L 
ANGE=ANG**(J-1 ) 

TF(ABS( ANGE) .LT.10.E-70) GO TO 7 
EM=EM+A ! J ) -ANGE 

E A 1 IS A FUNCTION WHICH DESCRIBES THE FEED VARIATIONS 
OVER THE PLANE OF THE LENS. 

E A 1=EM 

CER=CJ = AK*RAD*(COS ( ANG ) - ( COS (THETR)-l . 0 ) - S I N ( ANG ) * 

1 S IN ( THE TR ) * S I N ( A NG P ) ) 

THE SORT! L .O/COS ( ANG) ) TERM IS THE ONLY PART OF THIS 
EQUATION ASSOCIATED WITH THE LENS FUNCTION. 

THE REMAINDER IS FROM THE APPLICATION GF THE APERATURE 
FIELD METHOD TO THE SEMI-CIRCULAR APERATURE. TO MOD- 
IFY THIS PROGRAM TO PREDICT THE CONST ANT-K LENS FAR 
FIELD PATTERNS CHANGE THIS FUNCTION. THE CORRECT PHASE 
TERMS ARE ALREADY INCLUDED. 

C FCT4=C EXP !CER) *EA l-SQRT (COS { ANG ) ) 

RETURN 

END 
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